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ABSTRACT 



Context. Recent observations suggest that stellar magnetic activity may be influenced by the presence of a close-by giant planet. Specifically, 
chromospheric hot spots rotating in phase with the planet orbital motion have been observed during some seasons in a few stars harbouring hot 
Jupiters. The spot leads the subplanetary point by a typical amount of ~ 60° - 70°, with the extreme case of v And where the angle is ~ 170°. 
Aims. The interaction between the star and the planet is described considering the reconnection between the stellar coronal field and the 
magnetic field of the planet. Reconnection events produce energetic particles that moving along magnetic field lines impact onto the stellar 
chromosphere giving rise to a localized hot spot. 

Methods. A simple magnetohydrostatic model is introduced to describe the coronal magnetic field of the star connecting its surface to the 
orbiting planet. The field is assumed to be axisymmetric around the rotation axis of the star and its configuration is more general than a linear 
force-free field. 

Results. With a suitable choice of the free parameters, the model can explain the phase differences between the hot spots and the planets 
observed in HD 179949, v And, HD 189733, and t Bootis, as well as their visibility modulation on the orbital period and seasonal time 
scales. The possible presence of cool spots associated with the planets in r Boo and HD 192263 cannot be explained by the present model. 
However, we speculate about the possibility that reconnection events in the corona may influence subphotospheric dynamo action in those 
stars producing localized photospheric (and chromospheric) activity migrating in phase with their planets. 
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1. Introduction 
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About 300 extrasolar giant planets are presently knowrfl 
among which ~ 25 percent have a projected orbital semi-major 
axis lower than 0.1 AU. Such planets are expected to inter- 
act significantly with their host stars, not only through tides, 
but possi bly also with other mech anisms. Recent investiga- 
tions by lShkolnik etaD J2005U2008I) show that HD 179949 and 
v Andromedae have chromospheric hot spots that rotate with 
the orbital period of their inner planets. The spots are not lo- 
cated at the subplanetary point, but lead the planet by ~ 70° in 
the case of HD 179949 and by ~ 170° in the case of v And, 
respectively. They are quite persistent features with lifetimes 
of the order of 350 - 400 days, although they are not detected 
in all observing seasons. Two other stars, namely HD 189733 
and t Bootis, show some evidence of an excess of chromo- 
spheric variability, probably due to flaring, that is modulated 
with the orbital periods of their respective planets. Again, their 
active regions lead the subplanetary longitudes by ~ 70°. Wide- 
band optical photometry by MOST (the Microvariability and 



Oscillations of STar satellite; see, e.g., Wal ker et al . 2003) has 
given further support to the presence of an active region on r 
Boo leading the planet by ~ 70°. In 2004 it resembled a dark 
spot producing a light dip of abou t 0.001 mag, whil e in 2005 it 
varied between dark and bright (IWalker et al. 2008). 

Another intriguing observational evidence concerns 



HD 192263, a K dwarf for which ISantos et al.1 (120031) con 
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vincingly demonstrated the presence of a close-by planet with 
an orbital per iod of 24.35 days . Intermediate-band optical 
photometry bv lHenrv et al. (2002) showed that the star, during 
two observing seasons, presented a photospheric starspot that 
rotated with the orbital period of the planet for at least three 
orbital cycles. The spot was not located at the subplanetary 
longitude, but the planet was leading the spot by about 90°. 

In all the above cases, an interpretation based on tidal ef- 
fects can be ruled out because there is only one signature per 
orbit al cycle and not two. Following the original s ugges tions 
by ICuntz et al.1 d200dl) and iRubenstein & Schaeferl Q000), the 
chromospheric hot spot may be interpreted as an effect of the 
energy released by the reconnection between the stellar coronal 
magnetic field and the magnetic fie ld of the plan et that moves 
inside the star's Alfven radius (cf. lip et al.ll2004l) . To explain 
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the ph ase lag between the planet and the hot spot. lMcIvor et al 
(2006) determined the locations at which magnetic field lines 
that connect the planet to the star reach the stellar surface. At 
those locations, electrons accelerated at the reconnection site 
hit the chromosphere of the star producing the hot spot. Their 
model can account for the ~ 70° phase lag observed in HD 
179949 if the closed corona of the star extends out to the radius 
of the planet and the stellar large scale dipole field is signifi- 
cantly tilted with respect to the stellar rotation axis. However, 
the model cannot account for the ~ 170° phase lag observed in 
v And by any choice of the free parameters. 

A step towards a more realistic modell ing of the sta r -plane t 
magnetic interaction was made by ICranmer & Saarl (120071) . 
who extrapolated observed maps of the photospheric magnetic 
field of the Sun to compute the coronal field at different phases 
of activity cycle 22. The potential field source surface method 
was applied to derive the field assumed to become open at a 
distance of 2.5 solar radii from the centre of the Sun. By trac- 
ing magnetic flux tubes from the planet to the surface of the 
star, they were able to model the visibility of the hot spot along 
stellar rotation and the orbital period of the planet. The cyclic 
variations of the solar field accounts for the long-term varia- 
tions observed in the hot spot, while the inferred phase shifts 
between light curve maximum and planetary meridian passage 
range between -0.2 and 0.2. 



Preusse et alj ([2006) considered a different model in which 



the hot spot is due to the dissipation of Alfven waves produced 
by the planet moving inside the outer stellar corona. The key 
difference is that Alfven waves propagate along characteristics 
that do not coincide with the magnetic field lines and that can 
intersect the stellar surface virtually at any angle with respect 
to the direction of the planet, depending on the choice of the 
model free parameters. Therefore, their model is capable of ac- 
counting for both the phase lags observed in HD 179949 and 
v And. However, the chromospheric flux in a hot spot can be 
considered comparable to t hat of a typ ical solar active region, 



2 x 10 7 erg cm 2 s 1 (Priest 1982 ). which, in the case of 



i.e, 

isotropic emission, implies a very large Alfven wave flux at the 
source, located close to the planet. In order to have an energy 
flux at the source compatible wit h the estimates from the recon- 
nection model bv llp et alj (120041) . a highly collimated emission 
of Alfven waves is required, which is difficult to justify without 
a spe cific physical mech anism. Another difficulty of the model 
by iPreusse et a j] d2006h is that it cannot explain the observa- 
tions of t Boo and HD 192263 because Alfven waves cannot 
produce a cooling of the photosphere leading to the formation 
of dark spots. In view of such difficul ties, in the present w ork 
we move along the way indicated by iMcIvor et al. (2006), in- 
vestigating other geometries for the coronal field lines connect- 
ing the star with the planet and discuss the possible effects of 
coronal field reconnection on the dynamo action occurring in 
the convection zone of the star. 



2. Model of the stellar coronal field 

We focus on the interpretation of chromospheric spots moving 
synchronously with hot Jupiters. In Table[T| we list the relevant 
stellar and planetary parameters, i.e., from the first to the eighth 



column, respectively, the name of the star, its spectral type, ro- 
tation period P rot , radius R, semi-major axis of the planetary 
orbit (assumed to be circular) Aq, orbital period P or b, ratio of 
the Alfven velocity to the planet orbital velocity ma/moA, cal- 
culated for a magnetic field of 0.15 G and a coronal density 
at the orbital distance of the planet of 2 x 10 4 protons cm -3 ; 
and the references. We assume that magnetic reconnection at 
the intersection between stellar and planetary magnetospheres 
accelerates electrons to supra-thermal energies. A fraction of 
those high-energy particles travels inward along magnetic field 
lines and impacts onto the dense chromospheric layers of the 
star, producing a localized heating which gives rise to the hot 
spot. 

The position and the extension of the spot depend on the 
geometry of the stellar coronal field that conveys energetic par- 
ticles onto the chromosphere. We want to show that a coro- 
nal field geometry actually exists that can explain the observed 
phenomenology, although we shall not investigate how it can 
be produced and maintained. 

We assume that the orbital plane of the planet coincides 
with the stellar equatorial plane and that its orbit is circular. 
We consider a spherical polar reference frame with the origin 
at the centre of the star, the z axis along the stellar rotation axis 
and assume the radius R of the star as the unit of measure. We 
indicate the radial distance from the origin with r, the colatitude 
measured from the North pole with 9, and the azimuthal angle 
with (p. 

In the low-density environment of the stellar corona, the 
Alfven velocity is at least one order of m agnitude larger than 
the orbital velocity of the planet (cf. e.g., lip et al]l2004l and 
Table Q]), therefore we can assume that the corona is in a 
state of magnetostatic equilibrium. For the sake of simplic- 
ity, we consider a speci fic class o f magn etostatic equilibria, 
i.e., that investigated by iNeukirchl (1 19951) which includes as 
a particular case the linear force- free equilibria in sphe rical 
geometry previously described b y IChandrasekharl ([1956) and 
IChandrasekhar & Kendalll dl957l) . Generalizing a previous ap- 
proach by Low, Neukirch assumes that the current density in 
the corona has two components, one directed along the mag- 
netic field, as in force-free models, and the other flowing in a 
direction perpendicular to the local gravitational acceleration. 
Note that the gravitational field of the star dominates over the 
centrifugal force and the gravitational field of the planet over 
most of the stellar corona, because the star is rotating quite 
slowly (the rotation period of HD 1799 49 is about 7 days, 
while that of v And is about 12 days; cf. Shkolnik et alJl2005 . 



20081) and the mass of the planet is of the order of 10~ 3 stel- 
lar masses. The refore, the spherical equilibria computed by 
Neukirch! d 1995b under the assumption of a purely radial stellar 
gravitational field are good approximations for our case. 

The plasma pressure, density and temperature in the corona 
are not derived consistently by solving the energy equation, 
rather they are determined a posteriori to verify the magne- 
tohydrostatic balance and the equation of state. Nevertheless, 
in view of our ignorance of the plasma heating function, this 
approach is allowable to investigate analytically coronal field 
geometries of a wider class that linear force-free fields. 
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Star 


Sp. type 


Prat 


R 


A 




MA / Unrbi 
t\ i UtU 


Ref. fl 






(d) 


(Rn) 


(AU) 


(d) 






HD 179949 


F8V 


~ 7 


1.25 


0.045 


3.09 


14.6 


S08 


v And 


F7V 


~ 12 


1.25 


0.059 


4.62 


16.7 


S08 


t Boo 


F7IV 


3.2 


1.62 


0.046 


3.31 


15.3 


S08, D08 


HD 189733 


K1V 


11.7 


0.80 


0.031 


2.22 


15.2 


S08 


HD 192263 


K2V 


« 24 


0.78 


0.15 


24.35 


34.6 


S03 



" References: S08: Shkolnik et al. (2008); D08: Donati et al. (2008); S03: Santos et al. (2003). 



Non-axisymmetric coronal fields give rise to an additional 
modulation of the hot spot intensity with the rotation period 
of the star because magnetic field lines are rooted into the star. 
Therefore, we shall consider only axisymmetric configurations. 
Since the orbital radius of the planet is significantly greater than 
the stellar radius, the field component with the slowest decay 
as a function of the radial distance (i.e., the dipole-like com- 
ponent) is the one leading to the strongest interaction. For sim- 
plicity, we shall consider only that component assuming that it 
dominates over all the other components. This can be justified 
in the case of a solar-like a-Q. dynamo, as the axisymmetric 
dynamo mode with dipole-like symmetry with respect to the 
equator is the one with the strongest excitation rate and domi- 
nates the global field geometry outside the star. 

Under these hypotheses, we consider t he magnetic fiel d 
configuration indicated as case II in Sect. 4 o flNeukirchl (119951) . 
whose field components in the axisymmetric case (/ = 1, 
m = 0) are: 



R 2 

B r = 2B -rg(q)cos6, 

r z 

R 2 d 

B e = -B — — g(q) sin6, 
r dr 



B* = aB Q — g(q) sin0, 
r 



(1) 



where 2Bq is the surface magnetic field at the North pole of 
the star, a the helicity parameter of the corresponding linear 
force-free field, and the function g is defined as: 



g(q) = 



[^0-^-3/2(90) + coh/iiqo)] y/qo ' 



where bo and co are free constants, J -3/2 and 73/2 are Bessel 
functions of the first kind of order -3/2 and 3/2, respectively, 
q = \a\(r + a) and qo = \a\(R + a), with a a parameter measur- 
ing the deviation from the corresponding linear force-free field. 
The function g is a linear combination of the two independent 
solutions -yV + ^3/2(9) and V r + 1^3/2(9) = V r + a J -3/2(1) 
obtained by N eukirch for his case II (see Eqs. (40)-(43) in 
iNeukirchl [l995pl . When a — we recover the linear ax- 
isymmetric force-free solution o f Ichandrasekhai (1 19561) and 
Chandras ekhar & Kendall (1957) as given in, e.g. jFlver et al 



2 For the relationship between the Neuman n functions N n and the 
Bessel functions of the first kind J„, see, e.g., ISmimovl l fl964l) . 



( 2004). The physical meaning of the parameter a can be under- 
stood by noting that the transformation r + a — > r changes our 
field into a linear force-free field. Therefore, our field is a some- 
what radially "squeezed" version of a linear force-free field, 
with the additional confining forces being the gravity and the 
pressure gradient and a measuring the amount of radial com- 
pression. We introduce the absolute value of a in the definition 
of the radial variable q, so Eq. ([TJ is valid for both positive and 
negative a. Changing the sign of a is equivalent to change the 
sign of leaving B r and Bg u naltered, as can b e deduced by 
comparing Eqs. (18) and (24) of lNeukirchl (1 19951) . 

The magnetic field given by Eqs. (HJ does not extend to 
the infinity, but it is confined to the region between r - R 
and r — r^, where corresponds to the first zero of g (see 
Chandrasekhar 1956, for a discussion of the boundary condi- 
tions at r — rjj. Therefore, it can be used to model the inner, 
closed corona of the star, but not the outer region where the 
field lines become radial and the stellar wind is accelerated. 



Schriiver et alj d2003l) estimate that the extent of the closed 
corona is =s 2.5 for the Sun and increases remarkably for 
more active stars, e.g, as 19 for a star having a surface 
magnetic flux density or a coronal X-ray flux 10 times greater 
than the Sun. Since the stars we are considering rotate faster 
and generally have higher X-ray fluxes tha n the Sun (cf., e.g., 
Shkolnik et aill2005i [20081: ISaar et al.ll2007l) . we estimate that 
our model can be extended up to a radial distance of « 20 - 30 
before becoming inadequate (cf. Sect.O. 

To deri ve the path of th e magnetic field lines, we use the 
approach of iFlyer et al.1 (2004), introducing a flux function: 



(2) A(r,6) = B R 2 g(q) sin 2 0, 



(3) 



so that the magnetic field is given by: B = 
§7# + aA$|. Magnetic field lines lie over 
surfaces of constant A(r, 6), as can be deduced by noting that 
B ■ VA = 0. The variation of the azimuthal angle along a given 
field line can be derived from the equation: 



dr r sin Odfy 
B, Bj, 



(4) 



where dl = (dr, rd6, r sin 0d(f>) is the line element. Substituting 
the expressions of the magnetic field components, it gives: 



dq' 



cos 6(q') 



(5) 
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Making use of Eq. (01, and indicating by q e the value of q at 
which the field line reaches the equatorial plane, i.e., g(q e ) = 
sin 2 6>o, where 0q is the initial colatitude of the field line on the 
stellar surface, we find: 



<p(q) - </>(q ) 




g(q') 



g(q') - g(q e ) 



dq'. 



(6) 



Eq. (O is valid for a magnetic field line extending to the equa- 
torial plane because the integral exists all along the closed in- 
terval [go, q e ], as it is shown in Appendix lAl 

3. Applications 

The magnetic configuration connecting the chromospheric spot 
to the planet with the observed phase shift is not unique in 
the framework of our model, even in the force-free case. A 
complete exploration of the parameter space would be long 
and complex, so we limit ourselves to present some possible 
configurations obtained by fixing bo = -1.1, cq = 1.0, and 
a = 1.4, and varying a and 0q until the magnetic field lines 
reach the equatorial plane at the planet orbital radius r e with 
an azimuthal angle difference A0 e compatible with the obser- 
vations. Our choice of bo and cq is motivated by the need of 
having coronal field lines with an azimuthal angle A0 e suffi- 
cient to account for the large phase gap observed in v And. 
By making some experiments with different values of those pa- 
rameters, we selected a couple of values out of virtually infinite 
possible choices, only for illustrative purposes. In Table 12 for 
each star listed in the first column, we report the values of a, 
a, 6q, r e , A0 e , and in the columns from the second to the 
seventh, respectively. Note that corresponds to the first zero 
of g at q L = 3.609. For HD 179949 and v And we also fist the 
parameters for two force-free models. 

The best observed example of star-planet interaction is pro- 
vided by HD 179949, so we describe its models in some de- 
tail. A meridional section of its magnetosphere computed for 
a = -0.10 and a = 1.4 is shown in Fig.Q] The field lines at 
an initial colatitude 6q = 38.°85 reach the planet on the equa- 
torial plane at r e = 7.72 with A0 e comparable with the obser- 
vations - note that typical uncertainties on the observed phase 
lags are ±(15° - 20°). The corresponding force-free field with 
a = -0.10 gives G = 25.°0 and A<p e = 53° .7 which is too 
small to explain the observations. Therefore, instead of chang- 
ing bo or co, we decided to vary a and find an acceptable solu- 
tion with a = -0.12 that is also shown in Fig.Q] In this case, 
the field lines reaching the planet have an initial colatitude of 
#o = 26.°62. The angle A(p between the footpoint on the stellar 
surface and the point at distance r along a field line connecting 
to the planet is plotted vs. the radial distance in Fig. 12 At the or- 
bital radius of the planet this corresponds to an angle of -66. °4 
in the case of the field lines with a = 1.4, a — -0.10, giving a 
phase lag of ~ 0.18 between the planet and the hot spot, or to 
an angle of -73. °8 in the case of the force-free field lines with 
a = -0.12, giving a phase lag of ~ 0.205 . The latitudinal exten- 
sion of the hot spot on the surface of the star depends on the ra- 
dius of the planetary magnetosphere, i.e., the distance from the 
planet at which magnetic reconnection occurs. If we assume a 



magnetospheric radius R m = 4R p , where R p = 0.1 is the radius 
of the planet, then the radial boundaries of the magnetic recon- 
nection region are at 7.32 and 8.12, respectively. The colati- 
tudes on the stellar surface of the magnetic field lines reaching 
the equatorial plane at those distances are 39.° 17 and 38.°60, 
respectively. On the other hand, the longitudinal extension of 
the region is approximately 6°, as derived by the intersections 
of the surface of the star with the projections on the equatorial 
plane of the field lines tangent to the planet's magnetosphere. 
Therefore, the energy carried by high-energy electrons is fo- 
cussed onto a quite small region of the stellar chromosphere, 
producing a remarkable contrast with respect to the background 
emission. 

The hot spot on HD 179949 stays in view for about 0.5 
of an orbital cycle (see Fig. 6 in Shk olnik et al . 2008), which 
suggests that the inclination of the stellar rotation axis on the 
line of sight is quite high and the latitude of the spot quite low. 
Therefore, the magnetic configuration with a = 1 .4 is favoured 
with respect to the force-free configuration because it gives 
a spot at a lower latitude (~ 50° vs. ~ 65°). However, the 
non-force-free magnetic configuration is associated with sig- 
nificant perturbations of the thermal equilibrium of the stellar 
corona. Their relative amplitudes depend not only on the inten- 
sity and geometry of the field, but also on the background pres- 
sure and density distrib utions which are largely unknown (cf. 
Eqs. (28) and (30) of iNeukirchll 19951) . An illustrative example 
is shown in Fig. [3] where the relative variations of the tempera- 
ture, pressure and density in a meridional plane are plotted for a 
background isothermal corona with a field having the geometry 
computed for HD 179949. The background pressure and den- 
sity are proportional to exp[A(l/r- l/R)], with A = 10.4/?, and 
the ratio of the plasma pressure to the magnetic pressure at the 
base of the corona Bo ~ 5, while the unperturbed background 
temperature is 1.66 x 10 6 K. The largest temperature and den- 
sity perturbations are localized at high latitudes and quite close 
to the star. There the decrease of the plasma pressure provides a 
force that pushes the magnetic field lines toward the poles and 
the surface of the star, thus accounting for the squeezing char- 
acteristic of the considered non-force-free configurations with 
respect to the corresponding force-free models with a = 0, as 
discussed above. 

If the magnetic field of the planet is a dipole aligned with 
the planet's rotation axis and that axis is perpendicular to the 
orbital plane, then the planet's field is directed in the merid- 
ional plane close to the equator and reconnection occurs mainly 
with the stellar Bg component. If the star has a solar-like mag- 
netic cycle characterized by cyclic inversions of the sign of Bg, 
then reconnection delivers most of the energy when Bg is op- 
posite to the planet's field. This implies that the mean energy 
of the accelerated electrons is a function of the cycle phase of 
the star, with a significantly lower value when the two field 
components have the same sign. This may explain the on/off 
character of the star- planet magnetic interaction reported by 
IShkolnik et alJ (120081) . Recent observations of the evolution of 
the surface magnetic field in r Boo support the hypothesis that 
magnetic cycles in F-type dwarfs may be significantly shorter 
than the solar cycle (Do nati et aU 2008). 
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Models similar to that presented for HD 179949 can be 
computed for t Boo and HD 189733 for which the hot spots 
lead the respective planets by ~ 60° - 80°. On the other hand, 
in the case of v And, we need a higher value of a to account for 
the larger phase shift. In Fig. [4] we show meridional sections 
of the magnetosphere of the star computed for a model with 
a — 1.4 and a force-free model, respectively. Note the lower 
latitudes of the footpoints of the field lines connecting to the 
planet in comparison to the case of HD 179949. In Fig. [5] we 
plot the variation of the azimuthal angle along magnetic field 
lines for the models shown in Fig. [4] 

In the case of HD 192263, we may assume that chromo- 
spheric plages be associated with the photospheric spots lag- 
ging the planet by » 90°. The former can be explained by the 
present model by assuming a positive a, whereas in all the other 
cases the sign of a is negative (cf. Table [2]). Magnetic config- 
urations with opposite values of a have the same energy and 
opposite magnetic helicity, so there is no preference for a par- 
ticular sign of a from the point of view of the energy balance or 
field stability. Note that a linear force-free configuration with a 
given value of a 5 3 is stable provi ded that its m agnetic helic- 
ity exceeds a certain value (cf., e.g., Breg ell985h . 

We suggest that the sign of a may be related with the con- 
servation of angular momentum in radial plasma motions along 
field lines. In other words, if there is a mass flow along the field 
lines, it tends to bend the lines toward negative A(p (i.e., a < 0) 
when rotation is in the same sense of the orbital motion of the 
planet, or toward positive Acp (i.e., a > 0) when it is in the op- 
posite sense. The selection of the sign of a takes place close to 
the stellar surface since plasma motions cannot affect the field 
geometry in the low-/3 environment of the outer corona. This 
effect cannot be described in detail in the framework our mag- 
netohydrostatic approach, but it could be incorporated in more 
general models. 

It is important to note that in the case of HD 192263 the 
absolute value of a is 4 - 5 time smaller than in the other cases, 
because the relative orbital radius of the planet is correspond- 
ingly larger. As a consequence, the relative radius of the closed 
corona turns out to be 4 - 5 times larger than in the other cases. 
However, it can be significantly reduced by making a different 
choice for the free parameters bo and cq in the definition of the 
function g, in order to have its first zero at a value smaller than 
q L = 3.609. 



Magnetosphere of HD 179949 

n 1 1 1 p 1 




s/R 



Fig. 1. Projections of the magnetic field lines of the magneto- 
sphere of HD 179949 on a meridional plane. Configurations 
are displayed for Neukirch's models with two sets of param- 
eters, i.e., a = -0.10, a = 1.4 (solid lines) and a = -0.12 
and a = (dotted lines). The distance z from the equatorial 
plane and the distance s — r sin 8 from the rotation axis are 
measured in units of the stellar radius R. The small green circle 
indicates the planet assumed to have a radius of 0.1. The dot- 
dashed line indicates a field line with initial colatitude 38.°85 
connecting the surface of the star with the planet for the model 
with a = -0.10 and a — 1.4. The dashed line is the same for 
the force-free model {a = -0.12, a = 0.0) and has an initial 
colatitude of 25.°62. 

the toroidal field in a solar-like dynamo, i.e., the ao-effecjE 



4. Discussion 

The model presented in Sect. [2] can explain phase gaps up to 
~ 180°, whic h is not possib l e in th e framework of the previ- 
ous model by iMcIvor et al. (2006). Our model can accounts 
for chromospheric hot spots, but it cannot account for cool dark 
spots moving in phase with the planet, as suggested by MOST 
observations of t Boo in 2004 or in the case of HD 192263. 
Therefore, we speculate about the possibility that the magnetic 
reconnection in the stellar corona may perturb some process 
responsible for the dynamo action inside the star. Specifically, 
the effect leading to regeneration of the poloidal field from 



nection sites in the corona ( 


Blackman & Brandenburd 2003 ; 


Brandenburg & Subramanian 


2005). In Appendix |B~1 we pro- 



pose some preliminary ideas about a mechanism that may cou- 
ple coronal reconnection with the ffo-effect in the convective 
layers immediately below the surface. If a significant dynamo 
action takes place in thos e subsurface layers, as suggested by 



Brandenbur g ( 2005, 2007), this may give rise to a local amplifi- 



3 Note that the parameter qtd of the dynamo theory, considered in 
this Section, in Sect. [5] and in Appendix [B] is conceptually different 
from the parameter a characterizing force-free coronal field in Sects. [2] 
and[3] Therefore, we decided to mark it with the index D to avoid any 
confusion. 
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Fig. 3. Relative temperature (left panel), logarithm of the relative pressure (middle panel) and of the relative density (right panel) 
for an illustrative coronal model computed for HD 179949 (see text for details). 




Table 2. Model parameters for different cases of star-planet 
magnetic interaction 



Star 


a 


a 


e 




A0 e 


1 








(deg) 


(R) 


(deg) 


(R) 


HD 179949 


1.4 


-0.100 


38.85 


7.72 


66.4 


34.7 




0.0 


-0.120 


26.62 


7.72 


73.8 


30.1 


v And 


1.4 


-0.180 


48.90 


10.16 


173.6 


18.6 




0.0 


-0.208 


32.85 


10.16 


173.7 


17.3 


t Boo 


1.4 


-0.130 


44.50 


6.12 


72.7 


26.3 


HD 189733 


1.4 


-0.100 


38.48 


8.32 


77.9 


34.7 


HD 192263 


1.4 


0.025 


18.24 


41.19 


-106.7 


143.0 



Fig. 2. Variation of the azimuthal angle Atp along a magnetic 
field line starting at the surface of HD 179949 as a function of 
the radial distance in units of the stellar radius. The solid line 
is for the Neukirch's model with a = -0.10, a = 1.4, and an 
initial colatitude #o = 38°. 85; the dotted line is for a - -0.12, 
a = 0, and 6>o = 25°. 62. The vertical dot-dashed line indicates 
the orbital radius of the planet of HD 179949. 



cation of the magnetic field that will manifest itself with local- 
ized photospheric (and chromospheric) activity, including cool 
spots. Models with a longitude-de pend au effect have been ex- 
plored by, e.g jMoss et al. (2002) and show that the maximum 
magnetic field strength usually falls within ~ 15° - 20° from 
the maximum of the au effect. Differential rotation tends to 
wipe out non-axisymmetric components of the magnetic field, 
but active longitudes may survive at latitudes of minimal shear 
as discussed by, e.g., Big azzi & Ruzmaikin (2004). 

If this scenario is correct, we expect to find not only hot 
chromospheric spots but also cool photospheric spots at the 



longitudes where magnetic field lines connecting the planet 
with the star intersect the stellar surface. The phase lags ob- 
served in t Boo and HD 192263 may be explained in the frame- 
work of the coronal field model described in Sect. [3] 

In the case of the Sun, it is well known that new active re- 
gions have a marked preference to form close to existing active 
regions giving rise to complexes of activity w ith lifetimes up 
to 5 - 6 months (e.g., lHarvey & Zwaanlll993h . We conjecture 
that a possible cause of this increased probability of new mag- 
netic flux emergence may be a localized enhancement of the 
subphotospheric dynamo action promoted by helicity losses in 
the coronal part of an already existing active region. However, 
the correlation between the rate of coronal mass ejection and 
the rate of flux emergence in a complex of activity needs to be 
studied in detail to test this conjecture. 
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Magnetosphere of v Andromedae 




Fig. 4. Projections of the magnetic field lines of the magne- 
tosphere of v And on a meridional plane. Configurations are 
displayed for Neukirch's models with two sets of parameters, 
i.e., a = -0.18, a = 1.4 (solid lines) and a = -0.208 and 
a — (dotted lines). The distance z from the equatorial plane 
and the distance s = rsinfl from the rotation axis are mea- 
sured in units of the stellar radius R. The small green circle 
indicates the planet assumed to have a radius of 0.1. The dot- 
dashed line indicates a field line with initial colatitude 48.°90 
connecting the surface of the star with the planet for the model 
with a = -0.18 and a - 1.4. The dashed line is the same for 
the force-free model (a - -0.208, a = 0.0) and has an initial 
colatitude of 32.°85. 

5. Conclusions 

We extended the model proposed bv lMdvor et al. I d2006l) con- 
sidering a non-potential magnetic field configuration for the 
closed corona of a star accompanied by a hot Jupiter. We 
showed that a tilt of the axis of symmetry of the magnetic field 
with respect to the rotation axis of the star is not needed to 
reproduce the observed phase lags between the hot spots and 
the planets in HD 179949, v And, r Boo, and HD 189733. 
This agrees with mean-field dynamo models that predict that 
the large-scale stellar dipole-like field is aligned with the rota- 
tion axis when differential ro tation is significan t , as recently 



found in the case of t Boo (ICatala et al.l l2007t iDonati et alJ 
2008). The novelty of the present model is that it can easily 
account for a phase gap as large as that observed in v And, 




r/R 



Fig. 5. Variation of the azimuthal angle along a magnetic 
field line starting at the surface of v And as a function of the 
radial distance in units of the stellar radius. The solid line is for 
a Neukirch's model with a — -0.18, a = 1.4, and an initial 
colatitude 6 - 48°.90; the dotted line is for a = -0.208, a = 
0, and 6*o = 32°. 85. The vertical dot-dashed line indicates the 
orbital radius of the planet of v And. 



which is not possible with p revious approaches dMdvor et al 
200dlCranmer & Saarll2007l) . 



We also found that coronal linear force-free fields tend to 
give hot spots at quite high latitudes, making it difficult to 
explain the observed modulation of their visibi lity. Non-forc e 
free-equilibria, such as those computed by iNeukirchl ( 1995 ). 
give a significantly lower spot latitude for an appropriate choice 
of the parameter a, allowing us to fit the observations without 
assuming a high inclination for the stellar rotation axis and the 
orbital plane of the planet. 

The on/off nature of the star-planet interaction suggested by 
Shkolnik et al. (2008) can also be explained in the framework 
of our model as a con sequence of stellar activity cycles (see 



Cranmer & Saai 120071 for a model based on observed solar 



cycle variations of the magnetic field). 

We speculate about the role of magnetic helicity losses due 
to reconnection events in the stellar corona induced by the 
planet and suggest that this may give rise to a longitudinal de- 
pendent ffp-effect in the stellar dynamo process. Note that the 
mechanism we suggest may provide a physical basi s for the an 



effect perturbation model originally proposed by ICuntz et al 



(2000). Therefore, it is worthwhile to investigate dynamo mod- 
els with an an-effect varying in phase with the longitude of the 
planet to see whether they can reproduce the observations. This 
kind of models could provide an explanation for a modulation 
of the photospheric activity in phase with the planet that cannot 
be interpreted by means of the chromospheric heating model. 

Zeeman Doppler Imaging (ZDI) techniques are of great 
interest to map the magnetic fields in the photospheres of 
planet-harbouring stars, thus providing b oundary conditions 
for an extrapolatio n to their coronal fields ( Mo utou et alJl2007l : 
Donati et al.ll2008l) . Although non-potential configurations for 
assigned boundary conditions may not be unique, including 
ZDI information in our approach may give in principle the pos- 
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sibility of a self-consistent model of the stellar coronal field to 
study its interaction with the planet's magnetosphere (cf. e.g., 
Wiegelmann et al.ll2007l) . 
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Appendix A: On the integ rability of the function in 
Eq.© 

We prove that the integral appearing in: 



10) = \\ 



g(q') 



\ g(q') - gfae) 



dq'. 



(A.l) 



is finite, in spite of the denominator of the integrand vanishing 
at q = q e . Let us select a point q\ = q e — e, where e is a small 
positive number. Applying Lagrange's formula, we can write 
the denominator of the integrand as: g(q)-g(q e ) = g'(g)(q-q e ), 
where qi < q < q e ,q < £, < q e , and g' is the derivative of g. If g' 
has no zeros in the interval [qi,q e ], the function y/g(q)/g'(£;) is 
continuous and limited in the same interval. Indicating its upper 
bound with H, we have: 



~ f 

\ ■?(<?') -g(qe) ULl 2~X 

'♦!/• 

•Jci\ 



4>{q e ) - <p(q ) 
-dq' 



g(q') 



2J ?0 ^g(q')-giqc) 



dq' 



g(q') 



g(q') 



1 f 

1 r 

2 J qo \ g(q') - g(q e ) 

i C qi \ g(q') , , x „ i— 

2 1 im^) dq 



g(q') 



yg(q')-giqe) 

dq' 

yjq e - q' 



dq' 



qi- 



(A.2) 



This formula can be applied to estimate an upper bound for 
the error made when the numerical integration is truncated at 
a value q\ < q e . Specifically, in the case of the values listed in 
Table |2l we have fixed q\ by requiring that cos 6(q\) = 0.01, 
thus finding an error smaller than 3.° 6 for A</> e . 

Appendix B: Conjecture on a mechanism to 
perturb stellar dynamo action 

The value of the parameter au of the dynamo process depends, 
among others, on the mean valu e of the current helicity den- 
sity j ■ B (cf. e.g., Eq. (9.58) of Brandenburg & Subramanian 
120051) . Therefore, a perturbation of j ■ B in the subsurface layers 
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may produce a corresponding perturbation of an- The cur rent 
density in Neukirch's model (cf. Eq. (7) in lNeukirch1ll995 ) is: 



j = aB + VFx Vifr, 



(B.l) 



where ip is the gravitational pote ntial and F - F iVij/ ■ B,ifr) a 
suitable scalar function defined in lNeukirchldl995l) . Since VFx 
Vi^ = V x (FWif/) , taking the divergence of Eq. dB.lt we find 
that (B ■ V)a = 0, i.e., a is constant along a given magnetic field 
line. The current helicity at the surface of the star is given by: 



j ■ B = aB 2 + B ■ V X (FVifr). 



(B.2) 



If the perturbations of the coronal magnetic field due to recon- 
nection events do not directly affect the field at the surface, 
B 2 and B ■ V x (FVifr) are not modified by the motion of the 
planet. However, reconnection events induced by the planet 
will change the value of a in the corona and such changes will 
be conducted on the Alfven travel time to the stellar surface 
where they change j ■ B. Since the magnetic field at the sur- 
face is mainly radial, j ■ B ^ j n B n , where j n and B n are the 
components of / and B normal to the surface bounding the star. 
The perturbation of j n B n is then conducted into the subsurface 
layers because j n and B n are continuous across the bounding 
surface. 

According to this scenario, a localized subphotospheric 
perturbation of is produced, that moves in phase with the 
planet. A corresponding localized increase in the mean field 
gradient may promote a magnetic buoyancy instability that 
leads to the formation of an active region. Not e that magneti c 
buoyancy instability is a threshold process dAchesonl ll979). 
therefore a small increase of the field gradient may be sufficient 
to trigger the emergence of the magnetic field already amplified 
by the stellar dynamo. 

The active region will rotate with the orbital period of the 
planet, although the finite Alfven travel time across the corona 
and the finite time scale of magnetic flux emergence may intro- 
duce an additional phase lag between the planet and the spot. 
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